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There is mounting observational evidence that the expansion of our universe is 
undergoing an acceleration. A dark energy component has usually been invoked 
as the most feasible mechanism for the acceleration. However, it is desirable 



^ ■ to explore alternative possibilities motivated by particle physics before adopting 

such an untested entity. In this work, we focus our attention on an accelera- 
tion mechanism: one arising from gravitational leakage into extra dimensions. 
We confront this scenario with high- 2; type la supernovae compiled by Tonry et 
al. (2003) and recent measurements of the X-ray gas mass fractions in clusters 
of galaxies published by Allen et al. (2002,2003). A combination of the two 
databases gives at a 99% confidence level that flm = 0.291q;q2, ^rc = 0.2l;l^g;Qg, 
and Qk = -O.SGtoH, indicating a closed universe. We then constrain the model 
using the test of the turnaround redshift, Zq=o, at which the universe switches 
from deceleration to acceleration. We show that, in order to explain that accel- 
eration happened earlier than Zg=Q = 0.6 within the framework of gravitational 
leakage into extra dimensions, a low matter density, Q^n < 0.27, or a closed 
universe is necessary. 



Subject headings: cosmological parameters — cosmology: theory — distance 
scale — supernovae: general — X-ray: galaxies:clusters 
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1. Introduction 

The recent well known distance measurements of distant type la supernovae (SNela) 
suggest an accelerating universe at large scales (Riess et al. 1998, Perlmutter et al. 1999, 
Tonry et al. 2003, Harris et al. 2004, Knop et al. 2003, Riess et al. 2004). The cosmic 
acceleration has also been confirmed, independently of the SNela magnitude- redshift relation, 
by the observations of the cosmic microwave background anisotropies (WMAP: Bennett et 
al. 2003) and the large scale structure in the distribution of galaxies (SDSS: Tegmark 
et al. 2003a,b). It is well known that all known types of matter with positive pressure 
generate attractive forces and decelerate the expansion of the universe. Given this, a dark 
energy component with negative pressure was generally suggested to be the invisible fuel 
that drives the current acceleration of the universe. There are a huge number of candidates 
for the dark energy component in the literature, such as a cosmological constant A (Carroll 
et al. 1992; Krauss and Turner 1995; Ostriker and Steinhardt 1995; Chiba and Yoshii 
1999), a decaying vacuum energy density or a time varying A-term (Ozer and Taha 1987; 
Vishwakarma 2001), an evolving scalar field (referred to by some as quintessence: Ratra 
and Peebles 1988; Caldwell et al. 1998; Wang and Lovelace 2001; Weller and Albrech 2002; 
Gong 2002; Li et al. 2002a,b; Chen and Ratra 2003; Mukherjee et al. 2003; Gong 2004), the 
phantom energy, in which the sum of the pressure and energy density is negative (Caldwell 
2002; Dabrowski et al. 2003; Wang, Gong and Su 2004), the so-called "X-matter" (Turner 
and White 1997; Zhu 1998; Podariu and Ratra 2001; Zhu, Fujimoto and Tatsumi 2001; 
Alcaniz, Lima and Cunha 2003; Lima, Cunha and Alcaniz 2003; Feng, Wang and Zhang 
2004; Dai, Liang and Xu 2004), the Chaplygin gas (Kamenshchik et al. 2001; Bento et al. 
2002; Alam et al. 2003; Alcaniz, Jain and Dev 2003; Dev, Alcaniz and Jain 2003; Silva and 
Bertolami 2003; Makler et al. 2003), and the Cardassion model (Preese and Lewis 2002; Zhu 
and Fujimoto 2002, 2003; Sen and Sen 2003; Wang et al. 2003; Frith 2004; Gong and Duan 
2004a,b). 

However, the dark energy has so far no convincing direct laboratory evidence for its 
existence, so it is desirable to explore alternative possibilities motivated by particle physics 
before adopting such a component. In this respect the models that make use of the very 
ideas of branes and extra dimensions to obtain an accelerating universe are particularly 
interesting (Randall and Sundrum 1999a,b). Within the framework of these braneworld 
cosmologies, our observable universe is assumed to be a surface or a brane embedded in a 
higher dimensional bulk spacetime in which gravity could spread, and the bulk gravity sees its 
own curvature term on the brane which accelerates the universe without dark energy (Randall 
2002). Recently, based on the model of Dvah et al. (2000) of brane- induced gravity, Deffayet 
and coworkers (Deffayet 2001, Deffayet, Dvali and Gabadadze 2002) proposed a scenario 
in which the observed late time acceleration of the expansion of the universe is caused by 
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gravitational leakage into an extra dimension and the Priedmann equation is modified as 
follows 



where H is the Hubble parameter as a function of redshift z {Hq is its value at the present), 
flk, fire and flm represent the fractional contribution of curvature, the bulk-induced term and 
the matter (both baryonic and nonbaryonic) , respectively. Q^c is defined as Q^c = ^/^t^^Hq, 
where Tc is the crossover scale beyond which the gravitational force follows the 5-dimensional 
behavior. From a phenomenological standpoint, it is a testable scenario with the 
same number of parameters as a cosmological constant model, contrasting with models of 
quintessence that have an additional free function, the equation of state, to be determined 
(Deffayet et al. 2002). Such a possible mechanism for cosmic acceleration has triggered in- 
vestigations aiming to constrain this scenario using various cosmological observations, such 
as SNela (Avelino and Martins 2002; Deffayet, Dvali and Gabadadze 2002; Deffayet et al. 
2002; Dabrowski et al. 2004), angular size of compact radio sources (Alcaniz 2002), the 
age measurements of high- 2; objects (Alcaniz, Jain and Dev 2002), the optical gravitational 
lensing surveys (Jain et al. 2002) and the large scale structures (Multamaki et al. 2003). 
But the results arc disperse and somewhat controversial, with most of them claiming good 
agreement between data and the model while some of them ruling out gravitational leakage 
into an extra dimension as a feasible mechanism for cosmic acceleration. 

The purpose of this work is to quantitatively confront the scenario with the updated 
SNela sample compiled by Tonry et al. (2003) and to try to constrain the model parameters 
more accurately. It is shown that, although the two parameters, Qrc and Q^n, ai'c degenerate 
and there is a range on the parameter plane to be consistent with the SNela data, a closed 
universe is prefered by this scenario. As is well known, the measurement of the X-ray gas 
mass fraction in galaxy clusters is an efficient way to determine the matter density, fl^, and 
hence can be used for breaking the degeneracy between Q^c and Q^. When we combine the 
X-ray database pubhshed by Allen et al. (2002, 2003) for analyzing, we obtained a closed 
universe at a 99% confidence level, i.e., for the scenario of gravitational leakage into an extra 
dimension, a universe with curvature is favored by the combination of the two databases. We 
also analyze the turnaround redshift, -2^=0, at which the universe switches from deceleration 
to acceleration within the framework of the scenario. It is shown that, if the turnaround 
redshift happened earlier than Zq=Q = 0.6, only a low matter density, Q„i < 0.27, or a closed 
universe can explain this transition epoch. If, however, we consider the recent estimate by 
Riess et al. (2004), i.e., Zq=o = 0.46 ± 0.13, then a spatially flat scenario with = 0.3 (as 
suggested by clustering estimates) predicts Zq=o — 0.48, which is surprisingly close to the 
central value given by Riess et al. (2004). The paper is organized as follows. In the next 
section, we consider the observational constraints on the parameter space of the scenario 
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arising from the updated SNela sample compiled by Tonry et al. (2003), as well as the 
combination with the X-ray gas mass fractions in galaxy clusters published by Alien ct al. 
(2002, 2003). In section 3 we discuss the bounds on the model from the turnaround redshift, 
Zq=o. Finally, we present our conclusion and discussion in section 4. 

2. Constraints from SNela and galaxy cluster data 

Recently, Tonry et al. (2003) compiled a large database of SNela from the hterature 
and eight new SNcIa from the High-z Supernova Search Team. Since the techniques for data 
analysis vary between individual SNela samples, the authors have attempted to recompute 
the extinction estimates and the distance determination through the MLCS fitting (Riess et 
al. 1998), the Amis method of Phillips et al. (1999), the modified dml5 fitting (Germany 
2001) and the BATM method (Tonry 2003). Zero-point differences between each method 
were computed by comparing common SN measurements, and distances were placed on a 
Hubble flow zeropoint (dHo), and the median selected as the best distance estimate (for 
more details of this procedure, see Tonry et al. 2003; Barris et al. 2004). Tonry et al. 
(2003) present redshift and distance for 230 SNela, which includes many objects unsuitable 
for cosmological analysis, such as the SNela being heavily extinguished or nearby enough for 
velocity uncertainties to be a major problem. To determine cosmological parameters, the 
authors used a redshift cut of z > 0.01 and an extinction cut of Ay < 0.5 mag. The resulting 
sample of 172 SNela is illustrated on a residual Hubble Diagram with respect to an empty 
universe {Qm = 0, fire = 0) in Figure 1. We will use this sample to give an observational 
constraint on the model parameters, Qrc and Qm- 

For the ansatz (1), we are required to calculate as a function of the model parameters 
{flm, fire) and the Hubble constant Hq. Following Tonry et al. (2003), we added 500 km 
s~^ divided by the redshift in quadrature to the distance error given in their Table 15 
for calculating x^- order to concentrate solely on the density parameters, we need to 
marginalize over the Hubble constant Hq. Since Hq appears as a quadratic term in x^, it 
appears as a separable Gaussian factor in the probability to be marginalized over. Thus 
marginalizing over Hq is equivalent to evaluating at its minimum with respect to Hq 
(Barris et al. 2004). This procedure allows us to determine contours of constant probability 
density for the model parameters {flm, fire) corresponding to 68%, 95%, and 99% confidence 
levels, which is shown in Figure 2. The best fit happens at fl^n = 0.43 and fire = 0.26. As is 
shown in Figure 2, although there is a range on the parameter plane to be consistent with 
the SNela data, a closed universe is favored. Furthmore, the two density parameters, fire 
and flm, are highly degenerate, which is very similar to the degeneracy between fl\ and flm 
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found by Tonry et al. (2003). In order to determine Jl^c and Jl^ respectively, an independent 
measurement of Q^c or is needed. As shown below, the X-ray gas mass fraction data of 
galaxy clusters are appropriate for this purpose, because the data are only sensitive to flm 
(Allen et al. 2002, 2003). 

Since clusters of galaxies are the largest virialized systems in the universe, their matter 
content should provide a fair sample of the matter content of the universe as a whole, and a 
comparison of their gas mass fractions, /gas = Mgas/Mtot, as inferred from X-ray observations, 
with the cosmic baryon fraction can provide a direct constraint on the density parameter of 
the universe Q,n (White et. al. 1993). Moreover, assuming the gas mass fraction is constant 
in cosmic time, Sasaki (1996) shows that the /gas measurements of clusters of galaxies at 
different redshifts also provide a way to constrain other cosmological parameters describing 
the geometry of the universe. Recently, Allen et al. (2002; 2003) pubhshed the /gas profiles 
for the 10 relaxed clusters observed by the Chandra sateUite. Except for Abell 963, the 
/gas profiles of the other 9 clusters appear to have converged or be close to converging with 
r25005 the radius within which the mean mass density is 2500 times the critical density of 
the universe at the redshift of the cluster. The gas mass fraction values of these 9 clusters 
were shown in Figure 5 of Allen et al. (2003). This database can be used to break the 
degeneracy between fire and flm mentioned above, since it has been shown that the X-ray 
gas mass fraction is mostly sensitive to ilm no matter what the cosmological model is (Allen 
et al. 2002; Lima et al. 2003). The probabihty density over the model parameters, Qrc and 
Q^, for the 9 galaxy clusters is calculated using the method described in Allen et al. (2002). 
Following Allen et al. (2003), we include Gaussian priors on the bias factor, b — 0.93±0.05, a 
value appropriate for hot {kT > 5 KeV) clusters from the simulations of Bialek et al. (2001), 
on the Hubble constant, h = 0.72 ± 0.08, the final result from the Hubble Key Project by 
Freedman et al. (2001), and on fl^h^ = 0.0205 ± 0.0018 (O'Meara et al. 2001), from cosmic 
nucleosynthesis calculations constrained by the observed abundances of light elements at 
high redshifts. We then multiply the probability densities from the 172 SNela and the 9 
galaxy clusters, and obtain our final results on O^c and Qm, which are shown in Figure 3. 

Figure 3 illustrates the 68%, 95% and 99% confidence levels in the {Qm,^rc) plane using 
the red, green and yellow shaded areas, respectively. Our fits give at a 99% confidence level 
that Qrn = 0.291°:°^, Qrc = 0.21+°;°^, and Qk = -0.36+°J^. Although there is a range on 
the parameter plane being consistent with both the SNela and galaxy clusters data, and the 
resulting matter density is reasonable, a closed universe is obtained at a 99% confidence 
level, which is inconsistent with the result, fl^ — — 0.021q;o2, found by the WMAP (Bennett 
et al. 2003). Avelino and Martins (2002) analyzed the same model with the 92 SNela from 
Riess et al. (1998) and Perlmutter et al. (1999). Assuming a fiat universe, the authors 
obtained a very low matter density and claimed the model was disfavorable. In additional 
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to including new SNela data from Tonry et al. (2003), and combining the X-ray data of 
9 galaxy clusters, we relax the flat universe constraint in their analysis. We obtained a 
reasonable matter density, but a closed universe. In some sense, i.e., if we assume that 
our Universe is spatially fiat, as indicated by WMAP results, the accelerating scenario from 
gravitational leakage into extra dimension does not seem to be favored by observational 
data. However, two points should be emphasized here. First, that the same conclusion 
happens by performing a similar analysis with our current standard model, i.e., a ACDM 
universe. Second, that we have made heavy use of the X-ray gas mass fraction in clusters 
to determine the matter density. This kind of analysis depends on the assumption that /gas 
values should be invariant with redshift, which has been criticised by a minority of works in 
the field. For example, a recent comparison of distant clusters observed by XMM-Newton 
and Chandra satellites with available local cluster samples indicate a possible evolution of 
the M-T relation with redshift, i.e., the standard paradigm on cluster gas physics need to 
be revised (Vauclair et al. 2003). We should keep this point in mind when we make the 
conclusion that the gravitational leakage scenario is disfavored by the databases. 



3. Constraints from the turnciround redshift from deceleration to acceleration 

Since the scenario of gravitational leakage into extra dimensions is proposed as a pos- 
sible mechanism for the cosmic acceleration, the turnaround redshift from deceleration to 
acceleration is expected to provide an efficient way for verifying the model. It can be shown 
that the deceleration parameter as a function of redshift as well as the model parameters 
takes the form (Zhu and Fujimoto 2004) 

RR ^1 dlnS^ 

where E'^{z]Q.rc-i^m) — H^{z;^rc,^m)/ Hq. From Eq.(l), we could derive the turnaround 
redshift at which the universe switches from deceleration to acceleration, or in other words 
the redshift at which the deceleration parameter vanishes, which is as follows 

(l + .),.o = 2(^) (3) 



We have shown that Eq. (3) is generally valid no matter what the curvature of the universe 
is, though it was first obtained by Avelino and Martins (2002) for a fiat universe. According 
to Turner and Riess (2002), the value for the turnaround redshift lies in the la interval 
0.6 < Zq=o < 1.7. In Fig.4, the two dashed fines represent Zq=o — 0.6 and Zq=o — 1.7, 
respectively, while the hatched region at lower right corresponds to Zg=Q < 0, which means a 
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decelerating universe. The thick sohd hne is the flat universe. The vertical strip with cross- 
hatching is the matter density Q^n = 0.330 ±0.035 found by Turner (2002), and the vertical 
dot-dashed lines are Qm = 0.2, 0.4, a wider range. As is shown, in order to explain that 
cosmic acceleration started earlier than Zg=Q = 0.6, either a low matter density, Qm < 0.27, 
is needed on the assumption of a fiat universe, or a closed universe is necessary for a higher 
matter density. If, however, we consider the recent estimate by Riess et al. (2004), i.e., 
Zg=o — 0.46 ± 0.13, then a spatially flat scenario with fl^ = 0.3 (as suggested by clustering 
estimates) predicts Zq=o — 0.48, which is surprisingly close to the central value given by Riess 
et al. (2004). 



4. Conclusion and discussion 

The mounting observational evidences for an accelerating universe have stimulated re- 
newed interest for alternative cosmologies. Generally, a dark energy component with negative 
pressure is invoked to explain the SNela results and to reconcile the inflationary flatness pre- 
diction (fir = 1) with the dynamical estimates of the quantity of matter in the universe 
(Q^ ~ 0.3). In this paper we have focused our attention on another possible acceleration 
mechanism, one arising from gravitational leakage into extra dimensions. In order to be 
consistent with the current SNela and the X-ray clusters data, one would need a closed 
universe. 

Recently Lue et al. (2004) derived dynamical equations for spherical perturbations at 
subhorizon scales and computed the growth of large-scale structure in the framework of this 
scenario. A suppression of the growth of density and velocity perturbations was found, e.g., 
for n„i = 0.3, a perturbation of 5j = 3 x 10^'^ at Zi = 1000 collapse in the ACDM case 
at 2; ~ 0.66 when its linearly extrapolated density contrast is 6c = 1.689, while for the 
model being considered the collapse happens much later at z ^ 0.35 when its 5c — 1.656. 
Furthermore, the authors showed that this scenario for cosmic acceleration gave rise to a 
present day fluctuation power spectrum normalization (jg < 0.8 at a 2a level, lower than 
observed value (Lue et al. 2004). 

As is shown in Figure 2 of Deffayet, Dvah and Gabadadze (2002), on the assumption 
of a flat universe, luminosity distance for ACDM increases with redshift faster than that 
for the model being considered does (for the same Q^). Therefore it is natural that, if 
the ACDM model with {ilm = 0.3, il\ = 0.7, flk — 0) is consistent with the SNela data, 
the gravitational leakage model with {Qm — 0.3, fire = 0.1225, Qk — 0) will not be as the 
data are becoming enough to determine the cosmological parameters more precisely. While 
Deffayet et al. (2002) showed that the gravitational leakage scenario was consistent with the 
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54 SNela of the sample C from Perlmutter et al. (1999) - see also Alcaniz & Pires (2004) 

- Avelino and Martins (2002) claimed that this proposal was disfavored by the dataset of 
92 SNela from Riess et al (1998) and Perlmutter et al. (1999) [combining them via the 
procedure described in Wang (2000) and Wang & Garnavich (2001)]. We, however, think 
that only with a more general analysis, a joint investigation involving different classes of 
cosmological tests, it will be possible to delimit the flm — plane more precisely, as well 
as to test more properly the consistency of these senarios. Such an analysis will appear in a 
forthcoming communication (Alcaniz & Zhu 2004). 
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improved this work very much. Z.-H. Zhu acknowledges support from the National Natu- 
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(G2000077602). J. S. Alcaniz is partially supported by CNPq (305205/02-1) and CNPq 
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log z 

Fig. 1. — The 172 SNela data points, obtained by imposing constraints Ay < 0.5 and 
z > 0.01 on the 230 SNela sample of Tonry et al. (2003), are shown in a residual Hubble 
diagram with respect to an empty universe. The dashed and dot-dashed hues show (Q^, 
Q^J = (0.43, 0.26), our best fit, and (Q^, Qa) = (0.3, 0.7), the standard ACDM model, 
respectively. 
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Fig. 2. — Probability contours for fi^^ and fim in the model of gravitational leakage into an 
extra dimension, for the 172 SNela taken from Tonry et al. (2003) - see the text for a detailed 
description of the method. The 68%, 95% and 99% confidence levels in the - plane 
are shown in red, green, and yellow shaded areas, respectively. The cross-hatched region at 
the upper left represents the "no-big-bang" region, while the thick solid line corresponds to 
the flat universe. The best fit happens at Qm = 0.43 and = 0.26. 
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Fig. 3. — Probability contours over and fi^ for the combination of the 172 SNela taken 
from Tonry et ah (2003) and the 9 X-ray clusters from Allen et al. (2002, 2003). The 
68%, 95% and 99% confidence levels in the flr^ - plane are shown in red, green, and 
yellow shaded areas, respectively. The cross-hatched region at the upper left represents the 
"no-big-bang" region, while the thick solid line corresponds to the fiat universe. The best 
fit happens at flm = 0.29 and fire = 0.21, hence giving a closed universe with Qk = —0.36. 
However, the results depends on the X-ray gas mass fraction data from Allen et al. (2002, 
2003), in which the errorbars might be on the optimistic side. 
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Fig. 4. — Constraints on the parameter space (fire) ^m) of the model of gravitational leakage 
into an extra dimension from the turnaround redshift of acceleration. The hatched region 
at the lower right is the decelerating model, and the cross-hatched region at the upper left 
is the closed cosmological model without big bang. The right and left dashed lines are 
Zq=Q = 0.6, 1.7, respectively, while the thick solid line is the flat universe. Thus, in order to 
explain that acceleration happened earlier than Zq=o = 0.6, the gravitational leakage model 
needs a low matter density, Qm < 0.27, if the universe is flat. The vertical strip with cross- 
hatching corresponds to the matter density = 0.330 ± 0.035 found by Turner (2002), 
which clearly asks for a closed universe to explain Zq=Q > 0.6. For convenience, we also draw 
two dot-dashed lines for = 0.2, 0.4, for which there are some ranges to be compatible 
with a flat universe. We note that a matter density of < 0.27 is also permitted by the 



